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Proximity Effects near the Interface between d-wave Superconductors 

and Ferro/Antiferromagnets 
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We study the proximity effects near the interface between a d-wave superconductor (S) and 
a ferromagnet (F) or an antiferromagnet (AF). The S-side (F and AF-side) is described by the 
attractive (repulsive) Hubbard model, and the Bogoliubov de Gennes equation derived within 
the Hartree-Fock approximation is solved numerically. The superconducting order parameter 
and the magnetization (m) can coexist near the interface, and the spatial variation of m induces 
the spin-triplet p-wave component in both F and AF cases. The local density of states is also 
calculated and discussed. 
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Recently the proximity effect of unconventional su] 
conductors has been a subject of intensive study.B 
This is because the interface properties of these super- 
conductors can be quite different from those of conven- 
tional (s-wave) ones due to the nontrivial angular struc- 
ture of pair wave functions, so that their study is of par- 
tiucular interest. Now many unconventional supercon- 
ductors are known to exist, and the most famous exam- 
ples are the high-Tc cuprates in which_the (i3.2_j,2-wave 
superconducting (SC) state is realized.o' 

In this letter we study the proximity effect between a 
d-wave superconductor (S) and a ferromagnet (F) or an 
antiferromagnet (AF). We examine the possible coexis- 
tence of magnetism and superconductivity near the in- 
terface. In s-wave superconductors the coexistence is un- 
likely or limited because of the full gap in the excitation 
spectrum, while in the d-wave case it is not. We will show 
that the coupling of SC order parameters (OPs) with the 
magnetization induces a spin-triplet p-wave component. 
We also calculate the local density of states (LDOS) to 
show how the electronic structure changes across the in- 
terface. 

The system we consider is a two-dimensional d-wave 
superconducto(S)/magnet(M) bilayer. The direction 
perpendicular (parallel) to the interface is denoted as 
X (y), and we assume that the system is uniform along 
the y-direction except the two-sublattice spin structure 
in the AF case. We treat the Hubbard model on a square 
lattice within the Hartree-Fock (HF) approximation to 
describe the magnetism and superconductivity at zero 
temperature (T = 0). The crystal a-axis ia taken to be 
parallel to the x-direction, i.e., we consider only the [100] 
interface. The Hamiltonian for the two layers is given by 



Hl = -tL XI (4^cj, „ + ^-c-) + XI t^i^^T^a 

<i,j>cr i 

+ X ^^ [n^^n^i + ?^JT"a] > {L = M, S) 



where (?, j) and a denote the nearest- neighbor pairs and 
the spin index, respectively. Parameters t^, Ul and Vl 
are the transfer integral, the on-site interaction and the 
nearest-neighbor interactions, respectively, for the L{— 
M, S) side. The transmission of electrons at the interface 
is described by the following tight-binding Hamiltonian: 



Uj" — — tj" 



E 



{cla^rn,a + h.C.) 



(2) 



where I {m) denotes the surface sites of AI {S) layer, 
and then the total Hamiltonian of the system is H — 
Hm + Hs + Ht - fJ-J2ia Aa^io with [i bcing the chem- 
ical potential. We have examined the various values 
of tx/ti^usy The results are qualitatively similar, and 
the effect of proximity is reduced for smaller i^, as ex- 
pected. Thus the discussion is restricted to the case of 
It — tM = ts = t throughout this letter. 

The interaction terms are decoupled within the HF 
approximation 















(3) 

with Aij = {ci]Cji), and Aij and magnetization rrii = 
(jiii — nil)/ 2 ^^^ *^^ 0-f*^ ^o ^^ determined self- 
consistently. Along the y-direction we introduce two- 
sublattices A and B, and carry out the Fourier trans- 
formation: c^^ = ^fWlUk'^XA^Y''''' (" = A^?), 
where Ny (N^) is the number of sites along the y (x) 
direction. Now we define A" = 



A. 



(1)/ 






y i^^) = (cfrcf+si) and A);>{x,) = {cf^cf^^^). Then 
the mean-field Hamiltonian is written as (hereafter Xi is 
abbreviated as i) 



hmfa = XEE*lw^'^w*^w 



(4) 



<»j> 
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where 






F^ik) 

F,,ik) 
-W,,ik) 



W^j{k) 
F*,{k) 



Fj{k) 
-m,{k) 



ijl 



(fc) 



(5) 



*t(fc) ^ [4(fc)t,cf^(-fc),cf^(fc)t,cf,(-fc)]. (6) 



with 



W/,", = -t(<5,„+i + (5,,.-i) + (C/(<_,) - m)4 



l¥,,(fc) 



F,,(fc) 



— 2t cosfc (Jij 
-yAf,(<5,-.+i+<5,-,-i) 



1/ 



[e^'^■AW(^)-e-^■A(^)(^)]4 



(7) 

We diagonaUze the mean-field Hamiltonian by solvmg 
the fohowing Bogohubov de Gennes (BdG) equation:! 



y^ hi]{k)Ujn{k) = En{k)uin{k), 



(8) 



where En{k) and Ui„(fc) are the energy eigenvalue and 
the corresponding eigenfunction, respectively, for each 
k. The unitary transformation ^i(fc) = X^n ''^«"(^)rn(fc) 
diagonalizes the matrix TYmfa, and conversely the OPs 
Aij and rrii can be written in terms of En{k) and Ui„(fc). 
These constitute the self-consistency equations which 
will be solved numerically in the following^ 

In a uniform (bulk) case, the ground-state phase dia- 
gram of the repulsive Hubbard model ([/ > and V = 
within the HF approximation was examined by HirshJ3' 
There the ferro-, antiferro- and paramagnetic states are 
obtained depending on the value of U/t and the electron 
density. For an attractive case {V < and U = 0), 
the ground state is a spin-singlet SC state. Depending 
on the electron density a d^2_y2- (A^) or an extended 
s-wave (As)-wave SC state is stabilized, and the former 
is favored near half-filling. Here, A^ (A^) can be con- 
structed as a linear combination of A^'s in such a way 
that it changes its sign (is invariant) under 90° rotation 
in the basal plane. 

Now, let us study the M/S bilayer system. We im- 
pose the open (periodic) boundary condition for the 
X (y) direction, and the typical system size treated is 
N^ X Ny = 40 X 40 to 60 x 80 sites. We use t as a 
unit of energy (i.e., t = I). First we choose Um > 0, 
Vm = 0, Us — and T^s < so that the ferromag- 
netic and d-wave superconducting states are_iitabilized 
in M and S layers, respectively (F/S system).0' In Fig.l 
the spatial variations of OPs are shown. It is seen that 
the magnetization (m) and the SCOP (A) coexist near 
the interface, and the p^-wave component (Ap^;) as well 
as As appears. In general, the OP component different 
from that in the bulk can be induced due to the scatter- 
ing of Cooper pairs at an interface (or a surface faced to 
vacuum). The important point here is that we get the 
Pa^-wave state which is a spin-triplet (even parity) SC 
state. This state has pairing OP on the bonds along the 
x-direction. The spatial variations of OPs in the AF/S 
bilayer system are shown in Fig. 2. In this case the spin- 
triplet Pj,-wave component (Apj,) also appears. We will 



neglect the small A^ in the following discussions. 

The penetration of m into the S side causes the imbal- 
ance of the densities of spin-up and spin-down electrons. 
Thus electron pairs cannot be formed in singlet channels 
only, and then the spin-triplet component appears. It 
is seen that the M side is essentially unaffected by the 
proximity effect. This is because the Curie (or Neel) tem- 
perature is higher than T^ of superconductivity so that 
the typical lengths for the decay of the induced OPs are 
shorter in the M side. Then in the following we will 
mainly focus on the S side. 

We analyze the above results more precicely using the 
Ginzburg-Landau (GL) theory. (The GL theory is not 
quantitatively valid except near Tc, but it can give qual- 
itatively correct results.) The GL free energy in the S 
laver (in the continuum representation) can be written 
aE> 



^s = 



Kds 
Kdp 
K 



j—d,s,px,py 

{d,Ad){d.A 
Ad{{d.,m,)A; 



[a,\A,f + K,\dA,\'] 



■px 



{dyAa){dyAsy + c.c\ 
' {dym^)Apy}m^ + c.c\ 

'As{{d^m,)A;^ + {dym,)A;y}m^+c.c\). 

(9) 
This J-s is invariant under all symmetry operations for 
the square lattice and we have dropped higher order 
terms. All coefficients in J^s are positive definite ex- 



cept aj 

Vs' 



which are given at a temperature T as, aj 



1 



w 



(fc)tanh(Cfc/2T)/2a. for j ^d,s 



and aj = iVs/2) 



(Vs/N) Efe w^ik) tanh(6/2r)/6 
for j = Px,Py "Here N = N.j;Ny, ^k = -2t{coskx - 
cosky) — n, Wd{k) = coskx — cosky, Ws{k) — coskx - 



cos ky and w. 



'px{y) 



sink 



'x{y)- 



In the bulk superconduc- 



tor, the SCOP with the highest Tc occurs (Tc is given by 
a{Tc) — 0) and other components are usually suppressed 
due to the higher opier repulsive couphng to the bulk OP 
(not shown in jrg).E3^ Near half-filling Wd{k)'^ {ws{k)'^) is 
large (small) on the Fermi surface, and ad becomes neg- 
ative first. This is the reason we obtain the d-wave state 
near half-filling. (On the contrary, an extended s-wave 
state is favored away from half-filling.) 

In the M/S bilayer system, several OPs can coexist and 
they are not spatially uniform. In the F/S case, dxm ^ 0, 
dyTn = and A^ ^ 0, so that A^x is induced but A^y 
is not. On the contrary, in the AF/S system, finite dyvn 
induces Apy also. The opposite sign of A^^ and A^^ 
and their oscillating behavior in AF/S system can also 
be understood as being due to the staggered nature of 
m. The phases of the induced OPs relative to the bulk 
one (Ac;) are either or tt, since they are determined 
by the bilinear coupling terms in Ts- Namely, the state 
has the [dx^-y^ ±p2;)-syinmetry near the interface of the 
F/S system, and the (dx^-yi ipx ±pj,)-symmetry in the 
AF/S case. 

It has been proposed that antiferromagnetism and d- 
wave superconductivity jCaji^be treated in a unified way 
using SO{5) symmetry.tirEj' The present result seems 
consistent with this theory, since m penetrates deeper 
into the S side in the AF/S case than in the F/S case. 
On the contrary Apx decays faster in the former. The 
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reason for this is the following. The F/S system studied 
hereJs closer to half-filling compared with AF/S sys- 
tem.t3 Then, apx in J-g can be closer to zero (though 
still positive) in the F/S case because the factor w^^ is 
larger on the Fermi surface. The coherence length, which 
determines the decay of Aj , is given by S,j = ^JKj/oj^ 
and then the decay is slower in the F/S case. {Kj is not 
so sensitive to the electron density.) 

Recently, Honerkamp et al.liS studied the surface 
states of rf-wave superconductors. They found that 
near the [110] surface (faced to vacuum) the stag- 
gered magnetization can appear spontaneously, and that 
the p-wave as well as the extended s-wave compo- 
nent can occur. Their state breaks time-reversal sym- 
mptry spnntaripniisly arirl a surface Current is gener- 
ated.0BB@00MBe The reason the p-wave OP 
is obtained in the present case is, however, different from 
theirs. We consider a [100] surface and if it is faced 
to vacuum, the magnetization cannot be spontaneously 
generated, hence no Ap is obtained. In order to have a 
finite Ap for the case of a [100] surface, the proximity ef- 
fects from the magnets are necessary. Honerkamp et al. 
found the bound states near the Fermi level, which are 
different for the spin-up and spin-down components be- 
cause of local antiferromagnetism. If the [110]-oriented 
M/S bilayer is treated, the split bound states would also 
be expected. 

Next, jisffi show the results of the local density of states 
(LDOS).bl' The LDOS at site i{£ A) is given by 



Nf{i,u) 



NHi,oj) 



N. 2^ 2^ 



N, 



N, 



W4J-3,ri(fc)P'5(cj - En{k)) 



^^|u4z-2,„(fc)p,5(w + ^„(fc)) 



(10) 
(t,i being the spin indices) and NP{i,uj) (NP{i,uj)) is 
obtained by replacing (4i — 3) ((4i — 2)) in the above 
equation with (4i - 1), ((4z)). For the F/S case N^ = 
N^, and for the AF/S system N^ = N^^ {a =T, i)- In 
figures 3 and 4 (5 and 6), the results for the F/S (AF/S) 
system are shown. Deep inside the F (AF) layer, the 
LDOS for spin-up and the spin-down components is split 
as a result of magnetic orders (Fig. 3 (Fig. 5)). The spin- 
up and spin-down components gradually merge as the 
interface is crossed, and a T^-shaped DOS which is typical 
of a d-wave superonductor is formed. These changes of 
the LDOS as a function of the position appear smooth. 
In summary we studied the proximity effects near 
the interface between d-wave superconductors and 
ferro/antiferromagnets. We found the coexistence of the 
magnetic order and the superconductivity, which leads 
to the (d^2_j,2 ± px)-w&ve and (d^^-y^ zLp^ ±pj,)-wave 
surface states in the F/S and AF/S systems, respectively. 
The LDOS changes smoothly from magnetic to SC lay- 
ers. In this letter we did not take into account the effect 
of a vector potential that couples to the magnetization in 
the ferromagnetic state. This coupling may lead to more 
drastic effects, and is most suitably treated by field the- 
oretical methods. Such an effective theory can be con- 
structed based on the knowledge of the present solutions 
of BdG equations. This problem will be examined sepa- 



rately in the near future. 
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Fig. 1 Spatial variations of OPs for a F/S bilayer. Here 
Um = 10, Us = 0, Vm = 0, Vs = -2 and /i = -0.1. Note 
that all OPs are non-dimensional. 

Fig. 2 Spatial variations of OPs for an AF/S bilayer. 
Here Um = 4, C/g = 0, Vm = 0, Fs = -2 and ^ = 1. (a) 



Arf, As and m^(s). (b) Ap; 



A{B) 



Apy and A^. Indices A 



and B denote the sublattice. 

Fig. 3 LDOS for a F/S bilayer. The parameters are 

the same as in Fig.l and the finite width P = 0.08 is 
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introduced to each state, (a) x — —9a with a =|; (b) 
X — —9a with a —\,^ (c) x = 21a (a =t and | are almost 
degenerate). Here, a is the lattice constant and x = 
corresponds to the surface site of the M layer. 
Fig. 4 LDOS for a F/S bilayer. The parameters are the 
same as in Fig. 3 except x — a. (a) a =t, (b) a =J,. 
Fig. 5 LDOS for an AF/S bilayer. The parameters are 
the same as in Fig. 2 and the finite width F = 0.08 is 
introduced, (a) x = —9a with a =|; (b) x = —9a with 
a ~l, (c) X ~ 11a {a =t and I are almost degenerate). 
Fig. 6 LDOS for an AF/S bilayer. The parameters are 
the same as in Fig. 5 except x = a. (a) a =t, (b) a =[. 
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